Survival probability and order statistics of diffusion on disordered media 
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We investigate the first passage time tj^M to a given chemical or Euclidean distance of the first j 
of a set of A'^ ^ 1 independent random walkers all initially placed on a site of a disordered medium. 
To solve this order-statistics problem we assume that, for short times, the survival probability (the 
probability that a single random walker is not absorbed by a hyperspherical surface during some 
time interval) decays for disordered media in the same way as for Euclidean and some class of 
deterministic fractal lattices. This conjecture is checked by simulation on the incipient percolation 
aggregate embedded in two dimensions. Arbitrary moments of tj,N are expressed in terms of an 
asymptotic series in powers of 1/ In A which is formally identical to those found for Euclidean 
and (some class of) deterministic fractal lattices. The agreement of the asymptotic expressions 
with simulation results for the two-dimensional percolation aggregate is good when the boundary 
is defined in terms of the chemical distance. The agreement worsens slightly when the Euclidean 
distance is used. 

PACS numbers: 05.40.Fb, 66.30.Dn, 05.60.Cd 
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I. INTRODUCTION 



Diffusion in disordered media has been an area of inten- 
sive research during the last two decades |^ . Transport 
in non-crystalline, disordered materials cannot be ex- 
plained by the classical theories of diffusion since anoma- 
lous behaviors (relative to what happens in the Euclidean 
domain) are the rule here. In particular, the mean-square 
displacement of a random walker (r^(i)) is no longer pro- 
portional to the time t as occurs in uniform Euclidean 
systems of any dimension (Pick's Law), but the more 
general anomalous diffusion law (r^) ~ 2£'t^/'^™ holds 
for large times, r being the Euclidean distance from the 
position of the random walker at t = , the anoma- 
lous diffusion exponent, and D the diffusion constant. As 
the geometrical structure of real disordered media is very 
complex and varied, it is usually modeled by stochastic 
fractal lattices. Of these, the incipient percolation aggre- 
gate embedded in either two or three dimensions is the 
most widely used g |, |,|, |. 

Statistical problems related to a single random walker 
have traditionally been the subject of more intensive re- 
search than those corresponding to iV > 1 interacting 
or independent random walkers H, 0. Of course, prob- 
lems in which the walkers interact can not be analyzed 
in terms of the single-walker theory. However, there also 
exist other multiparticle problems that can not be ana- 
lyzed in terms of the single- walker theory even though the 
walkers are independent. These problems have begun to 
be the target of in-depth studies during the last decades 
|, |, 0, y, 0, 0, II, 0, |l|. Particular at- 
tention has been paid to (i) the evaluation of the average 
number of distinct sites visited up to time t hy N inde- 



pendent random walkers all starting from the same origin 
in both Euclidean and fractal lattices |8[ |ll|, p^ , 

and (ii) the description of the order statistic of the dif- 
fusion processes. This is the subject we address in this 
paper and can be stated as follows. A set of N indepen- 
dent random walkers all initially placed at a given site 
(the origin) of a medium start to diffuse at time t = 0. 
Eventually, at time ii,Ar(z) a random walker of this set 
reaches for the first time a site separated from the ori- 
gin by the distance z. Next, a second random walker 
reaches at time t2,N a site at the same distance z from 
the origin, and so on ||l^. Equivalently, we can under- 
stand tj,iq{z) as the time taken (escape time or lifetime) 
by the jth particle out of a set of TV 3> 1 to escape from 
a "spherical" region of radius z centered at the starting 
site of diffusion. Our goal in this paper is to calculate 
the escape-time moments (i^ jv(-^)) when TV 3> 1 random 
walkers diffuse in a disordered medium. 

This order-statistic problem was solved for the one- 
dimensional lattice in Ref. and for some class 
of deterministic fractal substrates in Ref. Por d- 
dimensional Euclidean lattices, the form of the first mo- 
ment of ti^N was guessed in Ref. ||l|| and checked using 
simulations results for d ~2 and d = 3. The full solution 
of this problem for Euclidean media and for arbitrary p 
and j has been obtained recently |jl|| . For all these me- 
dia, the pth moment of the time tj^^ spend by the jth 
random walker in reaching the Euclidean distance r was 
given in terms of an asymptotic series for large N of the 
form 
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where the coefficients T„m depend on p, j and the sub- 
strate. The question we want to answer is whether this 
is also true for disordered media. 

The prior knowledge of the short-time asymptotic ex- 
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pression of the survival probability r(z, i), which is de- 
fined as the probability that a random walker who starts 
at an origin site has not arrived at a spherical bound- 
ary of radius z in the time interval (0, i), was a key that 
allowed to find rigorously the full asymptotic approxi- 
mation (|l|) for (^j jv) ^'^'^ Euclidean and some determin- 
istic fractal media ^ in Refs. H, |l|, 0. (Here, and 
henceforth, we use the term "short time" or "short-time 
regime" to mean that 1, {z'^it)) being 

the mean square distance z travelled by a single random 
walker by time t.) However, the use of the asymptotic 
procedures of Refs. jl^, ^ for estimating (i^^jy) ^o^' 
disordered media is impeded by the fact that the value of 
the short-time survival probability is unknown for these 
substrates. In spite of this, in Ref. |Q it was conjectured 
that Eq. ([^) is valid for disordered media too |2l] . How- 
ever, Drager and Klafter using a scaHng approach, 
have shown recently that 



= o 



H(ln7V)i- 



(2) 



for ordered and disordered structures, which differs from 
(|I|) because, in general, d[^ ^ dyj. Here = dw/dmin is 
the chemical diffusion exponent defined by the relation 
{£) ~ const X t^/''"' , i{r) being the chemical or topological 
distance defined as the length of the shortest path con- 
necting two sites on a substrate that are separated by the 
Euclidean distance r, and c?min is the fractal dimension 
of the shortest path on the fractal: {({r)) ^ const x rf^^^ 
10, ^. Notice that dmin = 1 for Euclidean lattices and 
for the deterministic fractals discussed in Ref. so 
that, for these cases, Eqs. (0) and (||) agree. Unfortu- 
nately, the scaling approach of Drager and Klafter is not 
precise enough to lead to a fully correct asymptotic ex- 
pression of the form of Eq. (|^). This can be seen, for 
example, because the asymptotic series for (ti,Ar) given 
by these authors in Ref. jTs] does not agree with the rig- 
orous asymptotic series reported in Refs. |l^ for the 
one-dimensional lattice. 

In this paper we deal with the order statistic of 
the diffusion process in disordered media following the 
procedure previously carried out in Euclidean media 
and in (some class of) deterministic fractal 
Therefore, our first step must be to pro- 
pose a short-time asymptotic expression for the survival 
probability. The functional form we take coincides with 
that obtained for Euclidean and renormalizable fractal 
lattices. Next, we apply the asymptotic methods already 
used for the Euclidean and deterministic fractal cases 
in order to obtain asymptotic expansions for the order- 
statistics moments (i™^). 

The paper is organized as follows. In Sec. ^ a short- 
time asymptotic expression for the survival probability 
for disordered media is proposed. We check this conjec- 
ture in the two-dimensional incipient percolation aggre- 
gate and, on the way, estimate for this case the unknown 
parameters appearing in the proposed relationship. In 
Sec. HI we give the asymptotic expressions of the mo- 




ments of tj.A? for stochastic fractals assuming that the 
survival probability for the short-time regime is given 
by the expression proposed in Sec. ||. These theoretical 
results are compared with simulation data for the two- 
dimensional percolation aggregate. The paper ends with 
some conclusions and remarks. 



II. SHORT-TIME SURVIVAL PROBABILITY 



From now on, we will denote by z either the chemical 
distance i or the Euclidean distance r. As discussed in 
Sec. I we have to know the survival probability r(z, t), or 
equivalently the (boundary) mortality function h{z, t) — 
1 - T[z,t), for short times [i.e., for z / {z'^it))^/'^ » 1] in 
order to calculate the escape-time moments (i™iv)- This 
mortality function is the probability that a given random 
walker starting at the site z = will have been trapped 
by an absorbing boundary of radius z by time t. In this 
paper we will consider media for which this mortality 
function (averaged over all realizations of the lattice if 
the lattice is stochastic) grows for ^ = zjt^l'^'^ 1 as a 
stretched exponential with power-law corrective terms: 



/i(z, t) = /i(0 - exp her] { 1 + /iir" + ■ • ■ } 



(3) 



where A, /i, w, h\ are characteristic parameters of the lat- 
tice. The anomalous diffusion coefficient is replaced 
by dw if the Euclidean distance is used (z = r) and by 
d\^j\i z = i. There are good reasons to propose the func- 
tional form of Eq. (|3|) . To start with, this relation holds 
for Euclidean lattices. Then, the mortality function was 
obtained in Ref. jl4[ for some class of determinis- 
tic fractals by using a renormalization procedure which 
involved only boundaries containing the nearest neigh- 
bors of the origin after successive decimations, and the 
result is in agreement with Eq. (H) too. In these media 
z = r, fi ^ 1/2 and v = dw/{dw — 1). (Although the 
results of renormalization cannot be directly applied to 
any origin and any arbitrary boundary, there are again 
good reasons to suspect that the functional form is the 
same Q.) Also, a closely related quantity to h{z, t), the 
(site) mortality function ht{r) defined as the probability 
that a specific site r has been visited by a single random 
walker by time t, follows Eq. (^ for the two-dimensional 
percolation aggregate [^. Finally, one expects that the 
mortality function and the propagator P{z, t) decay for 
^ 3> 1 in the same way ]2^, and it is known that 
P(z, t) decays as cxp(— cf^) for stochastic fractals, where 
x; = d^/(d^ - 1) I, U ||. We will devote the remain- 
ing part of this section to verify that h{i, t) and /i(r, t) 
behave as conjectured in Eq. (^ for the two-dimensional 
percolation aggregate. 
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FIG. 1: Plot of ln[- In h{l,t)] versus ln| for the two- 
dimensional incipient percolation cluster. The radii of the ab- 
sorbing "circular" boundaries are i = 40 (triangles), £ = 100 
(circles) and £ — 160 (squares). The line represents the func- 
tion h{£, t) = il"*^* exp(-c|*) with i = 1, /i = -0.4, c = 0.9 
and V = dl/{di, - 1) = 1.714. 



A. Mortality function for the two-dimensional 
percolation aggregate: Chemical distance 

We start by estimating the mortality function h{£,t) 
when the "circular" absorbing boundary is placed at the 
chemical distance £ from the starting site. The numerical 
evaluation of this quantity is performed by the Chapman- 
Kolmogorov method (also called the exact enumeration 
method 0, |^). The circular boundary of absorbing traps 
is simulated by a set of special sites belonging to the clus- 
ter that absorbs all the probability density that enters 
them without giving back any probability to their neigh- 
bors. In the simulation we locate these boundaries at 
distances £ — 40, 100, 160, and perform three experiments 
of absorption for t = 0, 1, . . . , imax with t^ax = 1000 on 
every aggregate built, one for each boundary. The result- 
ing mortality function is averaged over 2000 realizations 
of the percolation aggregate, which are generated by the 
Leath method (|, |24| 

In Fig. we plot \n[—\nh{£,t)] versus In^ with ^ = 
£/t-^/'^w and = 2.4. [Hereafter we will put the sym- 
bol A ( '--^ ) over quantities corresponding to the chemical 
(Euclidean) distance.] The value of was taken from 
one of our previous works and is in agreement with 
the values reported by other authors ^ If the 
conjecture in Eq. (^ is right we should observe the lin- 
ear behavior ln[— \nh{£,t)]'^ Inc + -Oln^. Certainly the 
plots seem linear in Fig. ^ except for a portion in the 
range of large ^ where the curves deviate upwards. This 
is a finite size effect (already analyzed in the case of the 
two-dimensional Sierpinski gasket in Ref. |^) associated 
with the existence of a minimum arrival time correspond- 
ing to a random walker who travels "ballistically" along a 
chemical path from the origin to the absorbing boundary. 



1.0 1.5 2.0 2.5 

In I 

FIG. 2: Plot of ln[-ln/i(r,t)] versus ln| for r = 35 (trian- 
gles), r — 55 (circles), r = 75 (squares). The line represents 
the function h{r,t) = A^"^^ exp(-c|*) with A = l, fi = -1.8, 
c = 1.65 and v = d^/id^ — 1) = 1.56. 

which in turn implies a maximum available value of ^ in 
the simulations: ^max ~ £^^^^'^^. This value is around 
Imax - 2.15 for ^ = 40 and ^^ax - 2.96 for £ = 160. 
This apparently means that the reliable interval for nu- 
merical fitting is larger for ^ = 160 but we must also 
take into account that the minimum value of attained 

in the simulations is £/tll^x which is proportional to £ 
(in our simulations tmax = 1000). For this reason, in or- 
der to carry out the numerical fit, we have concatenated 
the simulation data for which the plots are almost lin- 
ear. The resuh h{£,t) = il"'^* exp(-cC*) with A ^ 1, 
c = 0.9, V = dL/idl - 1) ~ 1.714 and fi = -0.4 is also 
plotted in Fig.Il^. The agreement is excellent. 

B. Mortality function for the two-dimensional 
percolation aggregate: Euclidean distance 

We also simulated the mortality function h{r,t) in 
the percolation aggregate when the circular absorbing 
boundary has a Euclidean radius r. The analysis of these 
results parallels those of the previous subsection. In Fig. 
^ we have plotted the simulation results for the double 
logarithm of the mortality function ln[— lnft,(r, t)] versus 
In I (where ^ — r/t^^'^ ) for circular boundaries of radii 
r — 35, 55, 75, with t^na^ — 1000, and 2000 aggregates to 
perform the average. The anomalous diffusion coefficient 
for that time range was taken as = 2.8 in agreement 
with that obtained by us in a previous work |12[ . In the 
long time limit, a slightly greater value d^, = 2.87 has 
been found [Q, but as our simulations are restricted to 
t < 1000 it is more reasonable to consider the (effective) 
value dw — 2.8 which better represents the diffusive be- 
havior in this time range. 

In contrast with the results for the chemical boundary 
case, now we get a poorer collapse to a single line for 
the curves corresponding to different values of r. This 
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FIG. 3: Distribution of the mortality function with a circular 
trapping boundary, h{z,t), over 2000 realizations of the two- 
dimensional incipient percolation aggregate. We plot the his- 
togram J\f[ln h{z,t)] versus | In ^(z, t) | for t = 1000, z — r = 45 
(solid line) and z = £ = 80 (dashed line). The distribution 
is clearly wider for the Euclidean case than for the chemical 
case, although to a quite noticeably lesser extent than for the 
site mortality function [|l^ and the propagator j^, p5[ . 



we attribute to the fact that a given Euclidean distance 
r may correspond to many chemical distances £ > r 
depending on the site and the particular aggregate we 
are using in the simulation. The result is that, in con- 
trast with their chemical counterparts, the propagator, 
P{r, t) , the mortality function for a single trap site 
p2[ , and the mortality function for an absorbing spher- 
ical boundary, h{r,t) (see Fig. exhibit a broad dis- 
tribution over the different percolation aggregate real- 
izations. As a consequence, statistical precision requires 
more substrate averaging in the case of quantities defined 
in terms of Euclidean distance than quantities referred to 
the more natu ral c hemical distance. Anyway, proceed- 
ing as in Sec. IIB, we linked the simulation results for 
r = 35, 55, 75 where the plot is almost linear in order to 



perform an overall fit. Assuming the functional form of 
Eq. (H) we find that the data can be roughly described 



by h{r, t) ^ A^"'^^ exp(-£C* 
c = 1.65 and v = d^/^dw — 1) 



with A = I, jl 
~ 1.56 (see Fig. [ 



III. ORDER STATISTICS ON A STOCHASTIC 
FRACTAL LATTICE 

Now, with the functional form for h{z,t) found in the 
preceding section, we can proceed with the evaluation 
of (t™/v) by means of the asymptotic technique already 
developed for Euclidean and deterministic fractal sub- 
strates. Only the key steps in the calculation will be 
outlined. The interested reader is referred to previous 
references for the details ^ p^ . 

We will find it convenient in this section to write the 
mortality function as follows: 

h{t) EE h{z,t) - a(z)t^'^exp [-[to/tf] {l + hit'^ ...} 

(4) 

where a{z) = Az"'"', P = v/d^,, to{z) = c^/f^z< y hi = 
hiz~'" . The generating function of the mth moment of 
the jth passage time, UN,m{0 = Y,f=i{'f^J^N)C^~^ ^ can be 
written as lllSi 



1-C 



dtt"'-^ {[1 - h{t) + h{t)C\^ - c^} 



(5) 



Dropping the C, term we get 



poo 
77? / 

-— / exp {TV In [1-/1(0 (l-C)]} ■ 

^ ~ ^ Jo 

(6) 



The point is that UM.m{C} ^-i^d UN,miO have the same 
Taylor series expansion up to the term of order C 



N-l 



so 



that (tj^j^) can also be estimated through the evaluation 
of this pseudo-gcnerating function U at,™ (C) . Proceeding 
as in Refs. M, 14, |6|, one gets: 



UN,m{C) 



f 



1 



1 - C hi" A 



a (^ In In A — 7) 



2\nX 



(l + a) 



-2^(/i+(l + a)7)lnlnA + (l + a)/i2ln^lnA] + (^"yT^^ 



(7) 



where a = m/ (3 and 7 ~ 0.577215 is the Euler constant. 

Once the generating function is known, the escape 
times and their moments are calculated straightforwardly 
because (iJV) simply the coefficient of in the Tay- 



lor series expansion of UN,m{C)- Therefore, the mth mo- 
ment of the first passage time of the first out of ^ 1 
diffusing particles is equal to UN,m{0), i-e.. 
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ln"(Ao7V) 



a {fi In In XqN — 7) 
In XnN 



21n^(AoiV) 



(1 + a) ( y +7') +2/i7-2M^ 



-2/i(^ + (l+a)7)lnlnAoiV+(l + a)M^ln^lnAo7V] + O 



In'MnAoAf 
In^ XoN 



(8) 



with Ao = at^^ = Ac^ and a = m/ (3 = mv/df^ as before. 
Notice that the main term of this expression for m = 1 
agrees with the result given in [|l8) if we use v — df^/ (d^ — 
l). 

The calculation of (i"'jv) for > 1 is more in- 
volved. We begin with the identity ln"(l — z) = 
n! ^^jj— l)*5i(n)zYi!, where Si{n) are the Stirhng 
numbers of the first kind |^ . Using this relation, the 
expansions of 1/ In" A and similar terms in Eq. (|7|) in 
powers of 1/ In XqN can be found After some alge- 
bra we have 



where j = 2, 3, . 
A„(a) = H 



Nl 



ln"+i(Ao^) 



^ An(Q;) 

n 

n=l 



(9) 



. and 

a -I- 1 
In AoiV 



1 



■ 7 



O 



(n-1)! 
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In In XqN 

In^ XqN 



Finally, we will quote the main term of the asymptotic 
expression for the variance, ajj^ = (tJ^Ar) ~ {tj-N)"^, which 
is easily derived from Eqs. and ^) yielding 



2 



with 



'-0 



/32 (lnAoA^)2+2//3 



0-1 



dj = 



1+0 



In^ In Ao 
InAoA^ 



(11) 



(12) 



and j = 1,2,.... We use the convention that the sums 
are equal to zero when the upper limit is zero. It is 
clear that the main and first corrective terms of {t'j j^) 
are equal to those of (ij^Ar)^, so that only the difference 
between their second corrective terms contributes to the 
main term of the variance. For the sake of comparison 
with the simulation results it is more convenient to con- 
sider the quotient (tj n) /oj n whose expression is given 
by 



= PdjlnN 



1+0 



where di = 
(^2/6-5/4) 



(7rV6)-i/2, 

'1/2 



In-" \nN 
\nN ^ 

7rV6 - l)-i/2^ d: 



(13) 



Let us now compare the theoretical predictions in Eqs. 
(||) and ( p^ ) with simulation results in the two- 
dimensional percolation aggregate when an absorbing 
"circular" boundary is placed either at a given chemical 
distance z = £ or at a given Euclidean distance z = r. 



Absorbing boundary at a given chemical 
distance 



The first passage time of the first few random walkers 
out of a set of A^ = 2*, i = 2, 3, . . . , 16 to a "circular" 
boundary of chemical radius £ = 50 was simulated on 
2000 aggregates. In Fig. |^ we plot the scaled simulation 

results of [(ti.jv)/^'^™]"^ with 5 = l/(dt - 1), versus 
InA^ for £ = 50 and compare them with the theoretical 
predictions of Eq. (0) , namely. 
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5dl /i In In AT -7 
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In AT 
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where 5 = l/((i^ ~ 1)- Notice that 

pends linearly on In A^ if = dl^/{d^^ 
that this is indeed the case. 
From Eq. (^) one gets 



{{h,N)/^''-) ' de- 
- 1 ) . Figure ^ shows 



R 



{tl,N) 



t 



^ n 

\n=l 



V 



InAf. 



(15) 



so that R should be linear in In A^ and independent of 
the radius of the chemical boundary £. This is confirmed 
in Fig. H where this quotient is plotted versus In A^ for 
j = 2,3 and £ — 50,80. A good superposition of the 
simulation data for £ — 50 and = 80 is observed. The 
simulation results are in good agreement with the pre- 
diction of Eq. ([1) with V = dl/{dl - 1) and d^ = 2.4. 

Finally, in Fig. ^ we show the quotient between the 
average lifetime of the jth particle and its variance, 
{tj,N)/o'j,N, versus InA^ for the first, second and third 
random walkers to arrive at the chemical boundary with 
^ = 50. This is a remarkable figure that shows the 
crucial importance of the corrective terms on the order 
statistics formulae. Except for the first walker (j = 1) 
we find striking discrepancies between the zeroth-order 
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FIG. 4: Plot of T = [(ti,]v>/£'^™] versus InJV (iV = 2^ 2^ 

. . . , 2^'') for the two-dimensional incipient percolation aggre- 
gate where 5 = — 1). The circles are simulation results 
for £ = 50 and the lines are the zeroth-order (broken line) 
and first-order (solid line) asymptotic approximations with 
c = 0.9, V = 1.714 and /i = -0.4. 
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FIG. 6: Plot of kpj,N) joj^N versus In A'' for the two- 
dimensional incipient percolation aggregate with i — 50. The 
circles, triangles and squares denote the simulation results for 
j = 1, 2 and 3, respectively. The lines are the theoretical 
predictions with d^, = 2.4 and v = dl/(d^j, — 1) = 1.714. 
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FIG. 5: Plot of 7? = {t\,N)/{tj,N - ii.iv) versus InAf for the 
two-dimensional incipient percolation aggregate with j = 2 
and j — 3. The hollow [filled] symbols are the simulation 
results for ^ = 50 [£ = 80] which superpose closely except for 
large A'^. The lines correspond to the zeroth-order theoretical 
prediction \n{N)/[{di - l)Ei=i(l/")] with dl = 2.4. 



FIG. 7: Plot of In" In N/ In™ N versus the decimal logarithm 
LogioA'' for several pairs (n, m). These functions gauge the 
relative error of different finite-order asymptotic approxima- 
tions to the order-statistics quantities. 



(main term) of the asymptotic expression in Eq. dl^ ) 
and the simulation results. We attribute these discrep- 
ancies to the O [{In^ \n N) / \n corrective terms not 
considered in the zeroth-order asymptotic expression. 
Let us explain why. In Fig. ^ we plot the function 
In" InN/ In™ N for several values of n and m. One sees 
that In'^ In N/ In N cannot be neglected versus 1 even for 
values of N larger than the Avogadro's number. There- 
fore, it is not strange to find a poor agreement between 
the zeroth-order asymptotic expression and the simula- 
tion results. (In fact, what is truly unexpected is the 
good agreement of the zeroth-order theoretical prediction 
with simulation results for j = 1.) Notice that if we had 



calculated another corrective term in the asymptotic ex- 
pansion (^ , then the term neglected in Eq. (|l^) would be 
of order In* InN/ In^ N, which, as shown in Fig. ^ would 
likely be much smaller than the In^ In N/ In N corrections. 
Unfortunately, corrective terms of order In'' In N/ In N in- 
volve parameters such as hi [see Eq. (^)] which are very 
difficult to estimate relying only on simulation results. 
One also observes in Fig. that the relative errors of 
the zeroth- and first-order approximations to (tj,Ar) com- 
mitted by ignoring terms of the form In In N/ In TV and 
In^ In N / In^ N , respectively, will likely be small. This 
explains the good agreement between theory and simu- 
lation for lifetimes shown in Figs. and 0. 
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FIG. 8: Plot of r = [(ti,jv)/r''"'] versus IniV for the two- 
dimensional incipient percolation aggregate, where TV = 2^, 
2^, ... , 2^*^ and 5 = l/(d„ - 1). The filled [hollow] circles 
are the simulation results for r = 50 [r = 75]. The broken 
and solid lines are the zeroth-order and first-order asymptotic 
approximations, respectively, with c — 1.65, v = 1.56 and 
M = -1.85. 



B. Absorbing boundary at a given Euclidean 
distance 



We have also simulated the order-statistics of the ar- 
rival of N random walkers at a circular boundary of Eu- 
clidean radius r. The analysis of these results parallels 
those of the previous subsection. The numerical results 
for [ti,Ar/r'^™]-i/(d™-i) and {ti.N)/{tj,N - ti^) versus 
InA^ are plotted in Figs. || and|^, respectively. The radii 
of the circular boundary are r = 50 and r = 75, and 
we used 2000 aggregates to perform the substrate aver- 
age. The theoretical predictions of Eqs. ( p^ and ( p5| ) 
[replacing t by r, by d^,, 5 by l/{dyj — 1) and A by ^] 
lead to an agreement with numerical results that is not 
as good as that obtained when the chemical distance was 
used (see Figs. ^ and ||), although the linear behavior of 
T and R versus InA'^ and the superposition of the sim- 
ulation data for different values of r is still there. This 
slightly worse prediction is in accordance with the fact, 
discussed in Sec. II B, that the collapse of the numeri- 
cal mortality function for different distances r to a single 
curve is not so nearly perfect as when chemical distances 
are used. 



IV. CONCLUSIONS AND REMARKS 

In this paper we studied the order statistics of the sur- 
vival times (lifetimes or exit times) of TV ^ 1 independent 
random walkers that are put initially on a stochastic lat- 
tice at the center of a "spherical" absorbing boundary of 
radius z. We found the moments of the lifetime 

of the jth random walker of the set of N in terms of an 
asymptotic series that decays mildly in powers of 1/ In TV. 
These theoretical results were compared with numerical 



FIG. 9: The same as Fig. g but for circular boundaries with 
Euclidean radius r = 50 (filled symbols) and r = 75 (hollow 
symbols) . The lines correspond to the zeroth-order theoretical 
prediction ln(iV)/[(d„ - l)El=\(l/n)] with d» = 2.8. 



simulations for the two-dimensional incipient percolation 
aggregate with the absorbing boundary placed either at 
a given chemical distance z = ^ or at a given Euclidean 
distance z = r. The agreement found was reasonable 
and in accordance with the large size of the asymptotic 
corrective terms. 

The theoretical approach relies on the knowledge of 
the mortality function h{z, t) for short times. We as- 
sumed that the functional form of the short-time mor- 
tality function for a stochastic lattice averaged over all 
realizations coincides with that of the mortality function 
for Euclidean lattices and for the class of deterministic 
fractal lattices whose sites are isolated by their nearest 
neighbors. This is the main assumption of this paper and 
we showed that it is compatible with our simulation re- 
sults for the two-dimensional percolation aggregate with 
Euclidean and chemical absorbing circular boundaries. 
This allows us to describe the order statistics of the dif- 
fusion process on many kinds of substrates (namely, Eu- 
clidean media, and deterministic and stochactic fractals) 
with the same approach and formulae. 

A useful feature of the order statistics description is 
that it allows one to infer properties of the diffusive sys- 
tem (diffusion constant, number of diffusing particles, 
concentration of traps, effective dimension of the diffu- 
sive substrate, . . . ) from the analysis of the behavior of 
those particles that are trapped first. This could be an 
advantage when it is impractical or impossible to wait 
until the entire reaction is over. 

Finally, it should be noted that our simulation results 
for the stochastic fractal we have studied (the incipient 
percolation aggregate) arc insufficient to perform a com- 
pletely reliable numerical fit to the general form of the 
mortality function h{z, t) with its dominant exponential 
and subdominant power-law terms simultaneously. In- 
deed, even to find the true exponent value v of the dom- 
inant term for the propagator is not at all easy as past 
controversies about that true value reveal. Moreover, 



8 



the value of a hypothetical power-law subdominant term 
[equivalent to the term ^^^^ of Eq. (^)] is still a matter 
of discussion ||2^, Hence, one must expect that 

to find a completely reliable description of the short-time 
mortality function will be similarly difficult, and this task 
will therefore require more detailed simulations on larger 
lattices, with longer simulation times and, of course, av- 
eraging over many more realizations of the lattice than 
those we have used here. This is simply too much for 
a common workstation, but we think the effort on more 
powerful systems would be worthwhile: first, because of 
the interest in the percolation aggregate as a model of 
disordered media ||, and, second, because simulation 
of the mortality function (with circular boundary) is eas- 
ier than that of the propagator and, since the asymptotic 
coefficients and exponents of these two functions are re- 



lated, estimation of these parameters for the mortality 
function would help to solve the long-standing contro- 
versy on the exact asymptotic form of the propagator in 
fractal media. 
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